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Abstract 



Linearized Einstein gravity (with possibly nonzero cosmological constant) is quan- 
tized in the framework of algebraic quantum field theory by analogy with Dimock's 
treatment of electromagnetism [Rev. Math. Phys. 4 (1992) 223-233]. To achieve this, 
the classical theory is developed in a full, rigorous and systematic fashion, with par- 
\^q ■ ticular attention given to the circumstances under which the symplectic product is 

weakly non-degenerate and to the related question of whether the space of solutions 
is separated by the classical observables on which the quantum theory is modelled. 

o 

1 Introduction 

Low energy effects of quantum gravitation can usefully be modelled by the quantum field 
theory of linearized gravity, treated as a quantum field theory in a fixed background space- 
time. The background of de Sitter space is particularly interesting from the cosmological 
viewpoint, where linearized gravitons can induce fluctuations in the cosmic microwave 
background, for instance, see [5T1 Ch. 4 & 10] for a discussion of inflation and tensor fluc- 
tuations. For these reasons, the quantum field theory of linear gravitational perturbations 
has been studied extensively on several spacetimes, see for example [23] and references 
therein for the case of Robertson- Walker spacetimes. 

In the important de Sitter case, there has been a long-running controversy concerning 
the existence of a suitable de Sitter invariant vacuum state. This arises because natural 
constructions, for example, in the transverse-traceless and synchronous gauge associated 
with conformally flat coordinates on the Poincare patch [23] lead to a divergent expression 
for the two-point function as a result of an infrared divergence in the integral over modes. 
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Some parties to the controversy [33j nonetheless claim the existence of a de Sitter invariant 
two-point function whilst others [38J deny those claims. At issue is the validity of the use of 
Euclidean and analytic continuation methods, and the freedom to select gauge conditions 
or add gauge-fixing terms. 

For this reason it seems appropriate to have a framework for the quantization of linear 
gravitational perturbations on a general cosmological vacuum spacetime that will, when 
required, allow one to reduce down to a specific choice of spacetime and will also permit 
a rigorous investigation to be made into the issue of Hadamard states. This paper is the 
first step in that process. 

Of course the quantization of linear gravitational perturbations has previously been 
considered by numerous authors and it would not be feasible to list them all here, but, 
for our purposes [2] in particular stands out because of its consideration of Hadamard 
renormalization of the graviton stress tensor on cosmological vacuum spacetimes. However, 
this was achieved by introducing gauge-breaking terms and ghost fields. We prefer to avoid 
the introduction of such auxiliary fields and choices of gauge-breaking terms and apply a 
minimal approach along the same lines as Dimock's quantization of electromagnetism [16] 
(see also [17J). A brief outline of the method is as follows. Let (M,g) be a background 
spacetime that solves the cosmological vacuum Einstein equation. The classical phase space 
£P(M) of linearized gravity consists of solutions to the linearized Einstein equation modulo 
gauge-equivalence with a typical equivalence class denoted [7]. A classical observable is a 
function from &(M) to the set of complex numbers. For example, given any [7] e ^(M), 
we may smear any representative 7 (i.e., a solution to the linearized Einstein equation) 
against a compactly supported test-tensor /, to give an observable 



Ff(h}) = f la b f ab dv0l g , 



which is independent of the choice of representative provided that / satisfies the condition 
V a /( a b) = 0. Once a presymplectic form is given on &{M\ we may quantize to form a 
*-algebra of observables - essentially by Dirac quantization. This is a particular instance 
of the algebraic approach to quantum field theory in curved spacetimes (see, e.g., [IH1 
Ch. 4.5] for an introduction). The advantages of this procedure are: (a) it is manifestly 
independent of any choice of gauge, (b) it does not involve the addition of gauge-breaking 
terms or auxiliary fields in the action, (c) the method can be implemented in arbitrary 
globally hyperbolic cosmological vacuum background spacetimes, (d) it separates the con- 
struction of the algebra of observables from questions concerning the existence or otherwise 
of particular vacuum states which can then be addressed separately, and (e) it circumvents 
the known nonexistence of a Wightman theory of linearized gravity on Minkowski space 
allowing arbitrary smearings [4*6| [TO]. Nonetheless, there is a disadvantage, namely that 
we are not free to smear the metric perturbation against arbitrary tensors, which would 
be necessary in order to couple it to other fields. However, this quantization provides a 
rigorous framework that could be broadened in further work, for example, by exploiting the 
recent and rigorous approach to the Batalin-Vilkovisky formulation due to Fredenhagen 
and Rejzner [25J. Following the work reported here, Hack and Schenkel [30] have shown 



that our treatment of linearized gravity can be put into a broader framework of quantisa- 
tions of linear theories with gauge invariance. It would also be interesting to adapt these 
methods to other theories related to linearized gravity (T2J [13j IT4] . 

An important question is whether or not the class of observables identified above is 
sufficient to distinguish different equivalence classes of solutions. This is closely related 
to the issue of whether the pre-symplectic form is weakly non-degenerate, and thus sym- 
plectic. We will be able to resolve these questions positively at least in the case where 
the background spacetime has compact Cauchy surfaces, by adapting splitting results of 
Moncrief [39] derived in the Arnowitt-Deser-Misner (ADM) formalism (briefly reviewed in 
Appendix IBI). 

As mentioned above, the quantized theory leans heavily on the theory of classical lin- 
earized gravity and a substantial part of the present paper is devoted to a clear and general 
presentation of the theory. Much of this material is known, of course, but we have not found 
a unified and full treatment in the literature that would be sufficient for our purposesa 
We therefore hope that our presentation of the theory may be of independent interest and 
utility. To a large extent we emphasize a four-dimensional 'spacetime' viewpoint on the 
theory in contrast to the 'sliced' viewpoint of the ADM formulation. This is more natural 
for the formulation of the quantum theory and also removes any suspicion of dependence 
on particular slicings, coordinates or choices of (linearized) lapse and shift. However, as 
mentioned, we will draw on insights from the ADM formulation and make contact where 
necessary. 

The paper is structured as follows. In section [2] we begin with a brief introduction 
to linearized gravity and the gauge invariance of the theory. While our approach as a 
whole is independent of gauge choice, particular gauges are used for technical purposes 
to establish results on the full solution space and, to this end, we describe three gauges, 
de Donder, transverse-traceless and synchronous, and the circumstances under which they 
can be employed. This is largely standard, but for our discussion of the transverse-traceless 
gauge, for which we find that the topology of the Cauchy surface determines whether one 
may pass globally to the transverse-traceless gauge. In the synchronous case it is shown 
that a general perturbation is gauge-equivalent to a synchronous perturbation in a normal 
neighbourhood of any Cauchy surface (relative to the induced normal field). 

In section [3] we review existence and uniqueness of solutions to the linearized Einstein 
equation and provide the proofs for perturbations that are spacelike-compact, i.e., sup- 
ported in the union of the causal future and past of a compact set. Similar results are 
sketched in [22] for the case of compact Cauchy surfaces. We also prove various results 
concerning Green's operators, which will be required in section HJ The upshot of this 
discussion is that, modulo gauge equivalence, all spacelike-compact solutions to linearized 
gravity are given by the action of the advanced-minus-retarded solution, associated to the 
Lichnerowicz Laplacian P ab cd = V e V e S c a S d b — 2R c ab d , acting on smooth, compactly sup- 
ported, symmetric tensor fields / obeying V a {f a \, — \f c c g a b) = 0. One should note that 
although Lichnerowicz considers various Green's operators in [36J , his existence results are 
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purely local in nature, that is, they are valid in suitable open subsets of the manifold. 
Likewise, the book by Friedlander [26] also only considers local results. Our approach 
rests on the global theory of normally hyperbolic equations (developing the results just 
mentioned) that is described in the book of Bar, Ginoux and Pfaffle [6]. We make contact 
with Lichnerowicz's work by confirming that our propagator is identical to the one posited 
in [36]. He justified its form by using analogy with electromagnetism and the work of Fierz 
and Pauli [20] for the case that the background spacetime was Minkowski. In the present 
paper, it emerges (in any cosmological vacuum spacetime) as the Dirac quantization of 
Poisson brackets of the classical observables Ff. 

The final part is section HI which deals with construction of the phase space for lin- 
earized gravity and its subsequent quantization. As mentioned, the weak non-degeneracy 
of the symplectic product is established for compact Cauchy surfaces using results from 
the ADM formalism, in particular results of [39] providing various splittings of the space 
of initial data. Although these results are not established on general non-compact Cauchy 
surfaces, we conjecture that the non-degeneracy (as we formulate it) will hold for a large 
class of spacetimes with non-compact Cauchy surfaces as well. This is analogous to the 
situation in electromagnetism [16] where to prove non- degeneracy one appeals to the Hodge 
decomposition, again placing restrictions on the Cauchy surface. 

Once the phase space is constructed, the quantization then follows standard lines, either 
by Dirac quantization of the observables or a direct quantization of the symplectic space; 
we demonstrate that they are equivalent in section 14.31 We prove that the algebra satisfies 
the time-slice condition, i.e., it is generated by smeared fields whose smearing tensors are 
supported in any neighbourhood of a given Cauchy surface. We also briefly discuss the 
extent to which the construction respects the axioms of local covariance [TT] . 

Summarizing, our approach puts the quantum field theory of linearized gravity on a 
rigorous mathematical footing, which will permit future investigations into the states of the 
theory. In particular, it will enable precise statements to be made, using techniques from 
microlocal analysis, about Hadamard states (cf. [IT] for the case of electromagnetism). 
Elsewhere, we will investigate what can be said concerning the existence of a de Sitter 
invariant Hadamard state in the context of our framework. Note that the discussion of 
linearization instabilities would manifest itself in the properties of states [2] of the algebra 
of observables we construct and so we do not consider these issues here. 

We end this introduction with some preliminaries. We take, as in [311 08], the Riemann 
tensor to be defined by 

Kbc^d = (V a V b - V b V a K (1.1) 



and the Ricci scalar to be 



Rac - R a bc ■ (I- 2 ) 



Boldface type is used to indicate a tensor written with its indices suppressed. We also use 
w*° to denote the covariant form of a vector field w and v" to denote the contravariant form 
of a covector field v. A number of spaces of smooth scalar and tensor fields will be employed; 
all are taken to be complex- valued unless otherwise stated. As usual, C°°(M) denotes the 
space of smooth functions on M, while we write C°°(T^(M)) for the space of smooth rank 



(a, b) tensor fields on M. Various subscripts will denote restrictions on the support of such 
functions; C$°(M) and C^{T b a (M)) denote the compactly supported elements of C°°(M) 
and C°°(T b a (M)), while, following the notation of [6l p. 90], the subscript SC (spacelike- 
compact) attached to any of these spaces (e.g., C^ c {T b a (M))) denotes the subspace of 
tensor fields of the appropriate type whose support is contained within J(K) for some 
compact K C M. In particular this means that the intersection of the support with a 
Cauchy surface is compact. The subscript TC (timelike-compact) denotes tensor fields 
whose support lies between two Cauchy surfaces. Finally, S^M) (resp., S$(M)) denotes 
the symmetric elements of T®(M) (resp., Tq(M)), giving rise to corresponding spaces of 
symmetric tensor fields. 

2 Linearized gravity 

A spacetime is a pair (M,g) consisting of a four-dimensional, smooth, real, connected, 
Hausdorff, orientable manifold without boundary M together with a smooth Lorentzian 
metric g of signature ( — \- ++), with respect to which M is time-orientable; such a space- 
time is automatically paracompact by the theorem in the appendix of [29]. Each (back- 
ground) spacetime will be assumed to be globally hyperbolic, i.e., it admits no closed causal 
curves and for any two points p,q G M the set J + (p) fl J~ (q) is compact |8j. In addition, we 
assume that the metric solves the vacuum Einstein equations with cosmological constant 

G ab + Ag ab = 0, (2.1) 

where G ab = R a b — ^Rg a b as usual. As a consequence of fl2.ll) . the background Ricci tensor 
and Ricci scalar obey 

R ab = Ag ab , R = 4A, (2.2) 

which will often be used in what follows. 

Our aim is to quantize linearized perturbations of the Einstein equations. To linearize, 
we consider a smooth one-parameter family of solutions A h-> flf(A) to ( 12. ip with gr(0) = 
g. The linearized equation satisfied by the perturbation 7 := <jr(0) G C°°(T 2 (M)) on a 
cosmological background is 

Labfr) = 0, (2.3) 

where 

Lab{l) = --g a b(y c V d l{cd) - Dy - A7) - A7 (afe) - n 7(afe) 

-iv a V b7 + ^V c V (a7bc ), (2.4) 

7 = 7 a a is the trace of 7 and □ = g V Vf,. The linearized Einstein tensor L ab {^) vanishes 
for antisymmetric 7 and for symmetric 7 (our case of interest) reduces to 

Labh) = --g a biy c V d lcd - D 7 - A 7 ) - A 7a6 - -U lab - -V a V b 7 + V c V (a 7fe )c . (2.5) 



The linearized equation (12. 5 j) also follows from the Euler-Lagrange equations of the La- 
grangian (with 7 assumed to be symmetric) 

% = T abcde fV albc V dlef + S ahcd lablcdl (2.6) 

where T abcde f is given by 

rpabcdef _ _(q ad q bc q e S _|_ q a fq d ( b q c "i e -\- q d ( b q c ^ n ae 

- 9 ad 9 e{b 9 c)f ~ 9 a{e 9 f)d 9 bc ~ 9 d{b 9 c)a 9 ef ), (2.7) 

which is symmetric on interchange of b to c and on interchange of e to / as required by 
the symmetry of 7; T abcde f is also symmetric on interchange of abc to def. Finally 

nabcd „acbd , „bcad „ab ad fry q\ 

b =~^9 9 +^9 9 ~^9 9 (2.8) 

is symmetric on interchange of ab to cd and symmetric on interchange of a to b and 
interchange of c to d. The Lagrangian (I2.6P comes from the second order expansion of the 
Einstein-Hilbert action with cosmological constant 

S= I (R - 2A)dvol g , (2.9) 

where dvol g denotes the volume element on spacetime associated with g. Note that in the 
expansion, the linear term is a total divergence and the zeroth-order term is the Einstein- 
Hilbert action for the background. 

The covariant conjugate momentum 



1 A Q 

W bc = ^ = = 2T abcdefy ■ ^^ 

V^9 5 ^albc 



is given by 



uabc = _\ v a 1 bc + \ g bc v a 1 _\ g bc Vdl ad_\ g ac v b 1 _\ g ab v c 1+ \ v b 1 ac + \ v c 1 ab ^-Q 

La La La T: TT Lj £1 

and the Euler-Lagrange equations are thus 

V c IT cafe - 2S abcd ^ cd = L ab (j) = 0. (2.12) 

Note that the first equality holds for all smooth symmetric 7, not just for solutions to the 
linearized equations. 

The linearized Einstein equation (I2.3P is non-hyperbolic. To address existence and 
uniqueness properties of solutions we exploit gauge freedom to put the equation into a 
hyperbolic form. As is well known, the gauge transformations for linearized gravity take 
the form 

1 =7 + £w9, i-e-, 7o6 = lab + 2V (a u>&) (2.13) 



for w G C°°(T£(M)). Owing to the identity 

L ab (£ w g) = £ w (G ab + Ag ab ), (2.14) 

the linearized equation (12.31) is gauge invariant, in the sense that L ab (£ w g) = 0, if and only 
if the background is a cosmological vacuum solution (see [I5j Lem. 2.2]). This restricts our 
analysis to such spacetimes. 

In addition, to ensure well-posedness of the Cauchy problem and the existence of certain 
integrals, we restrict attention to spacelike-compact perturbations. With this in mind we 
introduce the notations 

&{M) = C? C (S° 2 (M)) (2.15) 

y{M) = { 7 G &{M) | L ab ( 7 ) = 0} (2.16) 

<f(M) = {£ w g | w g Cf c {Tl{M))} (2.17) 

for, respectively, the spacelike-compact symmetric tensors, the subspace obeying the lin- 
earized Einstein equation and the pure gauge solutions induced by spacelike-compact vector 
fields. 

In section H] it will be necessary to consider the larger class of spacelike-compact pure 
gauge solutions induced by arbitrary smooth vector fields, 

#(M) = {£ w g : w G C°°(T^(M))} n ST{M). (2.18) 

In fact, this coincides with @(M) whenever M has compact Cauchy surfaces; more gener- 
ally, the two sets differ only if there exist regions of the form M \ J(K), for K compact, 
supporting Killing vector fields. The example of Minkowski space shows that the latter 
condition, while necessary, is not sufficient: if £ w g G £f(M), then w G C°°(Tq(M)) satis- 
fies the Killing equation outside some set J(K) with K compact; without loss of generality, 
J(K) may be assumed to have connected complement, by expanding K if necessary. Be- 
cause Minkowski space is maximally symmetric, there exists a global Killing vector field 
£ such that supp(w; — £) C J(K); as £ w -^g = £ w g, we see that £ w g G £f(M), so 
£f (M) = £f (M) in this case. To see how £f (M) and 5f(M) can differ, consider Minkowski 
space with the causal future and past of the origin removed, which is still a globally hy- 
perbolic Einstein manifold M and inherits all the Killing vector fields of Minkowski space. 
Let r be the standard radial coordinate in the time-zero hyperplane and let K be the 
set of all points in the time- zero hyperplane with R < r < 2R, for some R > 0. Then 
M \ J{K) is disconnected, and neither component is relatively compact. Take any Killing 
vector field £ and let w = f(r)£, where / is constant outside [R,2R] with f(R) ^ f(2R). 
Then £ w g G £f (M) \ ^(M). 

The remainder of this section is devoted to identifying when particular gauge choices 
can be made, and is largely standard with the exception of our remarks on the transverse 
traceless gauge. We begin with the primary choice of gauge that will be used throughout 
this paper, the de Donder gauge, in which the linearized Einstein equations reduce to a 
hyperbolic form. 



2.1 de Donder gauge 

The de Donder gauge condition for a perturbation 7 is V a 7 afe = 0, where the standard 
trace-reverse operation (see e.g., |1EJ Ch. 7.5]) is defined by 7 afe := •y ab — \g a bl an d satisfies 
7 = —7 and 7 ab = jab- Note that for symmetric perturbations there is no ambiguity in 
writing V ■ 7 for V a, y ab , and we will often do so. 

To any rank (0,2) symmetric tensor, not necessarily a solution to fl2.3|) . there is a 
gauge equivalent rank (0, 2) symmetric tensor that satisfies the de Donder gauge condition. 
(Actually symmetry is not required in the proof.) Before stating this theorem we consider 
the following lemma concerning pure gauge perturbations. 

Lemma 2.1 For any w e C°°(Tq(M)), on a cosmological vacuum background spacetime, 
we have 

V-£Zg = (D + A)(wf. (2.19) 



Proof. By the definition of trace-reversal, (£ w g) a b = V ' a w b + V b w a — g ab V ' c w c . Taking 
the divergence of this gives 



V a (£ w g) ab = V a V a w b + V a V b w a - V b V c w c . (2.20) 

Using ( 12. 2 p we have V a V b w a = V&V a u> a + Aw b and hence the desired result. □ 

Theorem 2.2 The space ^(M) may be decomposed as 

ST{M) = ST dD {M) + Sf (M), (2.21) 

where 3? dD (M) = {76 ST{M) \ V a ^ ab = 0}. The intersection & dD (M) = ST AD {M) n 
^(M) is given by 



?dD 



(M) = {£ w g I w e C%> C (T*(M)), (□ + A)w = 0}. (2.22) 



Proof. Given 7 e 2T{M) let 7' = 7 + £ w g for an arbitrary w e Cffc(T$(M)). Taking 
the divergence of the trace-reversal of this and using Lemma 12.11 gives 

V a 7'a fe = V a 7 ab + (□ + A)i£%. (2.23) 

Therefore 7' G ^F d£) (M) if and only if w obeys 

(□ + A)w b = -V a 7?- (2.24) 

Applying standard results for hyperbolic equations, such as \6\ Thm 3.2.11] (generalized 
to permit non-compactly supported source terms as in Corollary 5 in |5, Ch. 3]), there is a 
unique solution to (j2.24p with initial data 10 |s and V n it>|£ that are smooth and compactly 
supported (we take them identically vanishing) on a smooth spacelike Cauchy surface S 



with future-pointing unit normal vector n. The solution is smooth and spacelike-compact c| 
Hence 7'— 7 G £f (M), which gives the splitting (j2.2ip . The final statement follows easily. □ 

Remark. From the proof, it is clear that the vector field w may be chosen to have support 
contained in the future (or past) of any given spacelike Cauchy surface. The space c S dD {M) 
specifies the residual gauge freedom within the de Donder class. 

Linearity of the equation of motion combined with the previous theorem gives: 

Corollary 2.3 The space 5?{M) decomposes as 

y{M) = y dD (M) + S?(M), (2.25) 

where S^ dD {M) = {76 3? dD (M) | L a \,{^) = 0} is the space of de Donder gauge solutions. 
Moreover, y dD (M) n &(M) = & dD (M). 

Next, we define a partial differential operator 

P ab cd :=n6 c a 5 d b -2R c ab d J (2.26) 

which is of the type considered by Lichnerowicz in [36] (see equation (10.4) of that refer- 
ence). We now note some important identities concerning this operator. 

Lemma 2.4 On a cosmological vacuum background spacetime, P commutes with trace 
reversal. In particular, P(t~) = if and only if P("y) = 0. 

Proof. We compute 



M - 1R c ah d )f cd = {U5 c a 5 d b - 2R\ d )f cd - -g ab {Uf + 2A/) 

= D(/ o6 - l -g ab f) - 2R c J(f cd - ±g cd f) (2.27) 

by using the identities g ab g ef R c e d f cd = -g a bR cd fcd = -Ag ab f and g cd R c ab f = -R ab f = 

—Ag ab f = g ab g e f R c e d f cd) which hold in cosmological vacuum spacetimes. □ 

Theorem 2.5 For any 7 G C 00 (5 , 2(Af )), on a cosmological vacuum background spacetime, 

2L afe ( 7 ) = -P ab cd l cd + (£^9) ab (2.28) 

or equivalently 

2LJ^~) = -P ab cd lab + {£(y.^9)ah- (2.29) 



b Here we use the fact that ~c has spacelike-compact support, which is not obviously the same as having 
support that has compact intersection with each Cauchy surface. After this paper was completed, and 
prompted by this concern, Sanders has shown that they are equivalent |44j . 



Proof. The Lie derivative term is 

7d\-7C 



'(V-tO 

jd-r-rc 



The P{"i) term is 



(£(v-y)WU = V a V c 7cb + V 6 V c 7c a - V a V 67 - g ab (y d V c lcd - -D 7 ). (2.30) 

One may show, using (J22D, that V a V c 7bc = V c V a7fec - A 7ab - R c ab d ^ c d and so 
(£{v.^9)ab = V c V a76c + V c V fe7ac - 2A 7a6 - 2i2 c a6 d 7«i - V a V b7 

-^(V d V c 7cd -^D 7 ). (2.31) 

- ^ Cd 7 crf = -0706 + ^D 7 + 2i? c a6 d 7cd + A<? afe7 . (2.32) 

Combine ( 12.311) and ( 12.32ft and then compare with ( 12.5ft . The second identity follows from 
the first by using Lemma 12.41 □ 

Therefore, for linearized gravity solutions on cosmological vacuum spacetimes that sat- 
isfy the de Donder condition V ■ 7 = 0, the equation of motion ( 12.4ft reduces to 

P ab c % d = D 7afe - 2R c ab d % d = 0, (2.33) 

or equivalently, by Lemma 12.41 

U lab - 2R c J lcd = 0. (2.34) 

For future reference, we note the following identity. 

Lemma 2.6 For any 7 e C°° (T® (M)) , on a cosmological vacuum background spacetime, 

V a (P ab cd lcd ) = (D + A)V a lab . (2.35) 

Proof. Expanding out the left-hand side of (12.351) gives 

V a (iV d 7cd ) = V a D 7ab - 2(V a R c ab d ) lcd - 2R c ab d W a lcd . (2.36) 

Substituting ( 12.2ft into the contracted Bianchi identity V a R bcd + ^bRcd — V c -Rm = gives 
V a R bcd a = and hence 

V a (P ab cd lcd ) = V a D lab - 2R c ab d V a lcd . (2.37) 

By using the Riemann tensor identity ( II. ip and the restriction to cosmological vacuum 
spacetimes ( 12. 2p one can show that 

V a U lab = UV a lab + AV a 7 a6 + 2R ad b c V dlac (2.38) 

and hence the result. □ 

It follows that if 7 is any solution to (I2.34p . its divergence, and the divergence of its 
trace reverse, obey 

(□ + K)(y a lab ) = (□ + A)(V a 7a6 ) = 0. (2.39) 

In addition, we see directly from ( 12.34ft that its trace obeys 

(D + 2A) 7 = 0. (2.40) 

In deriving equations (I2T33J) . ( l2T34p . (L2~39|) and (L2~40|) . the result (J2I2J was used. 
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2.2 Transverse-traceless gauge 

Many discussions of linearized gravity employ the transverse-traceless (TT) gauge, V a 7 a & = 
and 7 = 0, which is de Donder gauge with the additional constraint of vanishing trace. 
Inspecting the proof of Theorem 12.21 we see that 7 G 3T{M) can be put into the TT gauge 
if and only if the system 

(□ + A)u> a = 0, V a u> a = --7 (2.41) 

can be solved for w e C^ c (Tq(M)). It turns out that for vacuum spacetimes with a 
non-vanishing cosmological constant (e.g. de Sitter) this may be achieved at least when 7 
solves the linearized Einstein equation. 

Theorem 2.7 For cosmological vacuum spacetimes with A^O one may perform the fol- 
lowing decomposition of the space of spacelike- compact solutions: 

y{M) = y TT (M) + Sf (M). (2.42) 

As this is a departure from our main theme, the details are relegated to Appendix lAl 

By contrast, for vacuum spacetimes with a vanishing cosmological constant we find that 
there is a cohomological obstruction to the solution of ( 12.41 ft : it is possible if and only if 
the trace of the solution 7 obeys 



V n 7rfS = (2.43) 

on some (and hence all) Cauchy surfaces. To some extent this problem arises because we 
only consider spacelike-compact perturbations. If this restriction is dropped and £ is non- 
compact, then the TT gauge can also be achieved for A = 0, as is the case in Minkowski 
space. 

2.3 Synchronous gauge 

The synchronous gauge is defined relative to a spacelike Cauchy surface S by the condition 
n a, ~Yab = 0, where n is the future-pointing unit normal vector to S. Here, we describe how 
one can make a gauge transformation to put any solution into the synchronous gauge near 
a Cauchy surface. The result is similar to jU Lem. 1.1], which shows that this can be done 
on a Cauchy surface; our result is therefore slightly more general and our proof treats the 
solution of various equations arising in detail. 

Before we state the theorem we recall some definitions. Given a submanifold S C M, 
the normal exponential map exp -1 is the restriction of the exponential map to the normal 
bundle of S; therefore all the geodesies arising from this map will meet S orthogonally. A 
normal neighbourhood of S is a neighbourhood of 5* that is diffeomorphic under exp x to a 
connected neighbourhood of the zero section in the normal bundle of S. In particular, on 
any normal neighbourhood & of £ there is a unique future-pointing, geodesic, hypersurface 
orthogonal, unit vector field that we call the normal field of E in iff. 
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Theorem 2.8 Let E be a smooth spacelike Cauchy surface with future-pointing unit normal 
vector n. Let G be any open normal neighbourhood of E, whose closure is contained in 
another normal neighbourhood o/E. Then 

&{M) = ^ ch (M) + #(M), (2.44) 

where ^ y ^ c (M) = {76 S"(M) \ n a ^ ab = on G} and n is the normal field o/E in G. 
In particular, n a ^ a b = holds on E. 

Remarks. E has normal neighbourhoods by j42j Prop. 7.26]. Given any such normal neigh- 
bourhood, we may restrict to a smaller normal neighbourhood whose closure is contained 
in the original. Therefore the existence of G in the hypothesis is not restrictive. 

Proof of Theorem 12.81 Let 7 G £?(M) be arbitrary. Then the condition 7 + £ w g G 
S^ V C (M) amounts to the equations 

V a W = -\h a h b lab (2.45) 

(V^||) 6 - w\\ a V b fi a = -n a lab - V b W - -n a h c lac h b , (2.46) 

where Wo = n a w a and w\\ — w + W^h . 

The first step is to obtain a solution Wq to ( 12.45ft on & . This can be achieved as 
follows. Through each point p G E we have a unit speed normal geodesic X p : I — > @ with 
G / C K and A p (0) = p. Equation (I2.45P can be integrated along \ to give a solution 

1 /"* 
(W o X p )(t) = -- / {n a n h iab o \ p ){s)ds (2.47) 

and we define the scalar function Wq at any q G G by W / o(q') := (W^o ° A p )(t g ), where p g 
and t q are uniquely determined by q = X Pq (t q ). This satisfies ( 12. 45 ft by definition. It is 
smooth on O because t q and p q vary smoothly with q under the normal exponential map, 
which is a diffeomorphism on 0, and n a h b, y a b is smooth by assumption. 

Now, to obtain wn on G we solve ( I2.46P locally within a neighbourhood of each geodesic 
and then patch together the results with a partition of unity. For each q G E, let jV q C E be 
an open normal neighbourhood of q. Hence on jV q we have well-defined normal coordinates 
x % (i — 1, 2, 3) based at q and associated basis vector fields ej. Now, for each g G E let ^ g 
be the open set of points in G connected to E by geodesies emanating normally from jV q . 
The sets ^# g , for q G E, form an open cover for G because it is a normal neighbourhood. 
Within each ^l q we can introduce Gaussian normal coordinates given by: the proper time 
t along the geodesies, with t = on E, and the normal coordinates x % mentioned above. 
In these coordinates, (12.461) becomes 

d(w\\)j(t,x) j dW (t } x) 

j t 2r J . (t,x)(w||),(t,x) = - 7w —. . (2.48) 
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This system can be solved using standard results (see, e.g., [47J Sec. 1.6]) to give w\\ on 
^ q , where we have also used that (itf||)o = in these coordinates. This process is repeated 
on each ^ q for all q G E. 

As E is an embedded submanifold of M, it will also be second-countable and Hausdorff. 
Therefore by [501 Thm 1.11], the open cover {^V q \ q G E} of E by normal neighbourhoods 
will admit a countable partition of unity {xx | A G /} subordinate to the cover with supp xx 
compact for each A G /. Hence, for each A G / there exists a q G E such that supp Xx C. jV q . 
To obtain a suitable partition of unity xx on the ^# g 's, we solve VnXx — with Xx\s = Xx 
by integrating along integral curves of n as before. 

Therefore the w\\ that we seek on O is given by 



W\\ 



J^Xa^J, (2.49) 



where each wb is the solution to ( I2.48P on the set jtft q that contains xx- Observe that 
( I2.49P satisfies (I2.46P on & by the properties of xx- 

In conjunction with Wq this will give the w on & to transform to the synchronous 
gauge. We now examine the support properties of w. Outside supp 7, ( I2.45P reduces to 
VnWo = and so Wq = constant along each normal geodesic emanating from E, as long 
as the geodesic does not enter supp 7. Choosing Wq\t, = yields Wq = on every geodesic 
that does not intersect supp 7; hence Wo\e is spacelike-compact. Using this means that 
outside supp 7 equation ( I2.46P reduces to (VnW\\)b — w\\ a Vbn a = 0; in Gaussian normal 
coordinates, the right-hand side of (I2.48P vanishes. Thus with 1011 |s = 0, the solution w\u 
vanishes in every ^# g that does not intersect supp 7, so w\ff is spacelike-compact. Outside 
& we let w smoothly decay to zero. Hence w may be chosen to be compactly supported 
and therefore £ w g G @(M). □ 

3 Existence and Uniqueness of solutions 

We will now prove existence and uniqueness up to gauge of solutions to the non-hyperbolic 
equation (j2.4p . and show how Green's operators of the de Donder hyperbolic equation 
(j2.34p are related to linearized gravity solutions. 

3.1 Existence 

Let E be a smooth spacelike Cauchy surface with unit future-pointing normal vector n. 
The Cauchy data mart] Data E : C^iT^M)) -> C oo ((T 2 (M)| s )©C oo ((T 2 (M)| s ) is defined 
by 

Data s (7):=(7| s ,V n 7|s)- (3.1) 

Our aim is to solve the equation L ah {-y) = subject to given Datas(7). As is well known, 
this cannot be achieved for arbitrary data, because the components -^(t^Ie do not 



c The notation T^M)]^ means the restriction of rank (0, 2) tensor fields on M to the surface S. 
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involve second time-derivatives of 7 and are therefore completely determined by the initial 
data. Put more formally, there is a constraint map C s : C oo ((T 2 (M)| E )©C oo ((T 2 °(M)| s ) -»> 
C°°( T 0(M)| s ), so that 

C s ( 7 |s,V Tl 7|s):=n a L ab (7)| s (3.2) 

and the Cauchy data must be restricted to the kernel of C . The precise form of the linear 
map C E will not be needed here. We observe that gauge invariance of L ab entails that 
C s o Datas is also gauge invariant. 

Theorem 3.1 Let E be a smooth spacelike Cauchy surface with future-pointing unit normal 
vector n. For any initial data C £ £ Cq° (S^M)^) satisfying the initial value constraint 
C S (C,£) = there exists a solution 7 £ 3?{M) to {23} such that Data s (7) = (£,£). 

Proof. The proof is broken into two steps. First, we make a gauge transformation that 
puts the initial data into the de Donder gauge on E. Second, as the initial value problem 
for (j2.34p is well-posed, we obtain a solution with the transformed data, which obeys 
the de Donder condition globally; it is here that the initial value constraint is vital. By 
Theorem 12.51 the solution will satisfy the linearized Einstein equation ( 12. 4p and so by 
undoing the original gauge transformation we obtain a solution with the original Cauchy 
data. 

Following the structure just set out, we begin by constructing a smooth x £ &(M) 
such that Datas (x) = (CO- This is accomplished as follows: In a normal neighbourhood 
of E, we use parallel transport along geodesies normal to E to obtain £ £ 3F{M) with 
Datas(C) = (CO)- Taking any extension of $, in 3F{M\ we form \ = £ + s£, where s 
is determined uniquely at each point p of the normal neighbourhood by the requirement 
p = exp q sn\ q for some q £ E. By extending smoothly, we obtain x £ 2T{M) with the 
required properties. 

Using the splitting of Theorem l2.2l there exists w £ C^ c (Tq(M)) such that 7 = x+£w9 
obeys the de Donder condition V a 7 a6 = and 

(%(, - 2V( a Wfe))|s = Xab|s = Cab 

n c V c (j ab - 2V( a w b) )| s = n c V c Xab|s = £a&- (3.3) 

Passing to the second step, let 7 £ £T{M) be the (unique) solution to the hyperbolic 
equation Dt^ — 2R c ab d ^f cd = with initial data Datas(7) — Datas(7). The existence and 
uniqueness of 7 follows from [HI Thm 3.2.11]. 

The key point is now to show that 7 obeys the de Donder condition globally. As 7 obeys 
(I2.34p then by Lemma I2.6[ V a 7 a6 obeys the hyperbolic equation (I2.39J) ; it also vanishes on 
E because V a 7 ao |s = V a 7 afe |s = 0. The following identity is now required. 

Lemma 3.2 On a cosmological vacuum background spacetime, for any solution 7 to (j2.34p . 

it holds that n c V c (V a 7 ao )|s = 2L ab (j)n a \ s . 
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Proof. Combining the hypothesis with Lemma I2T41 and using Theorem 12.51 gives 2L ab {^f) = 
(£(y.-s\t9)ab- Contracting with n and expanding the right-hand side gives 

2L a6 ( 7 K = n a V a V c ^ cb + n a V b V c ^ ca - n b V d V c ^ cd . (3.4) 

The metric may be written in terms of the normal vector n and a projection operator 
(see [3U Ch. 2.7]), so that g ab = —n a n b + q ab . This allows one to split a vector into its 
components normal and tangential to E. Using that V a 7 afe |s = and n a n a = — 1 one finds 
that 

2L afe (7K| s = n c V c (V^ a6 )| s . (3.5) 

As stated immediately prior to Theorem 13. 1[ the constraints are gauge invariant and 
so we have the following chain of equalities 

MtKIs = C s (Data s ( 7 )) = C s (Data s (7)) = C s (Data s (x)) = (3.6) 

and thus n c V c (V a 7~ afe )| s = 0. 

Accordingly, we have shown that V a 7 afe obeys a hyperbolic equation with vanishing 
initial data on E; it therefore vanishes globally in M by [6], Cor. 3.2.4]. Thus, 7 solves 
(j2.34p and satisfies the de Donder condition, so 7 G y dD (M). By undoing the original 
gauge transformation, we obtain a solution 7 = "f — £ w g to the linearized Einstein equation 
f 12 .51) . which obeys Datas(7) — (C>£) by virtue of ( 13. 3ft . 

3.2 Uniqueness 

Given any initial data satisfying the constraints then by Theorem 13 . 1 1 there exists a solution 
to the linearized Einstein equation. However, as the next theorem shows, this solution is 
only unique up to gauge equivalence. 

Theorem 3.3 Suppose 7,7' G S^iM) with Data^(7) = Data?: (7') on some spacelike 
Cauchy surface E. Then 7 = 7' + £ w 9, also written as 7 ~ 7'. If, additionally, 7,7' G 
S^ dD (M) then gauge equivalence is replaced by equality. 

Proof. Let £ = 7 — 7' which satisfies Datas(^) = 0. By Theorem 12.21 we may write 
7 = £ + £ w g where 7 G y dD (M) and w G Cf^T^M)) obeys (□ + k)w b = -V a £ a & . 
Choose initial data, io|s = and Vntu|s = 0. Therefore Datas(7) is given by 

7afe|s = (V a w b + V 6 w a )| s = (3.7) 

™ c V c 7Je = n c V c (V a w b + V b w a )| s . (3.8) 

In fact n c V c 7 a b|s also vanishes, as we will now show. Firstly, by the choice of data we know 
that V a Wb\-£ — and therefore any derivative of this taken tangentially to E will vanish. 
Next, using the Riemann tensor identity (11.11) and the conditions on w at E, we have 

n c V c V a w b \x = V a (n c V c w b )\j: - (V a n c )V c u> 6 | s + n c R cab d w d \^ (3.9) 

= V a (n c V c w b )| s = -n a n d V d {n c V c w b )\x. (3.10) 
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As Datas(^) = 0, we have (n + A)iu|s = — V-£|s = 0. Expanding, using g a b = —n a nb + q a b, 

n a V a (n c V c w b )\v = q ac V a V c w% + {n a V a rf)V c w% + Kw b \^ = 0. (3.11) 

Therefore n c V c V a Wb|s = and hence by (13. 5J) , n c V c ^ a b\i: vanishes. As 7 satisfies the hy- 
perbolic equation (12.341) with vanishing Cauchy data, it vanishes globally by [6, Cor. 3.2.4]. 
Thus £ ~ and hence 7 ~ 7'. 

If both of the solutions 7 and 7' are de Donder with the same initial data then they 
must coincide because they solve the hyperbolic equation (I2.34[) with identical data. □ 

3.3 Green's operators 

Any de Donder solution satisfies the two equivalent hyperbolic wave equations (I2.33P and 
(I2.34p . By the results of \6\ Ch. 1.5], the differential operator P (defined in equation (I2.26P ) 
is a normally hyperbolic operator and therefore, by [61 Cor. 3.4.3], admits unique ad- 
vanced (-) and retarded (+) Green's operators E ± : C^(T$(M)) ->■ C°°(T$(M)), whose 
action on test tensors we can write as 



(E ± ab cd f c , d ,)(x)= / E^o'ixrffMWvolW, (3.12) 

J M 

where / G C °°(T 2 (M)). The operators E ± satisfy P{E ± f) = E ± P(f) = f for all 
/ G (^(Tg (M)), and have the support properties snpp(E ± f) C J^supp/); moreover 
^ = E ± f is the unique solution to -P(7 =t ) = / with support that is compact to the 
past (+)/future (— ) (for the uniqueness of the Cauchy problem see [HI Thm 3.2.11]). The 
advanced-minus-retarded solution operator is defined to be E := E~ — E + . Analogous 
properties hold for all normally hyperbolic operators. Any solution to (I2.34[) may be written 
in terms of E as the next lemma shows. 

Lemma 3.4 Any 7 G ^(M) solving P(7) = may be written as 7 = Ef with f G 
CrOSSC^- 
Proof. See P Thm 3.4.7]. □ 

The preceding lemma concerns general solutions to the hyperbolic equation (I2.34[) . how- 
ever, as we are looking for solutions to linearized gravity we are interested more particularly 
in those solutions satisfying the de Donder condition, V-7 = 0. Before showing when such 
a solution will be de Donder, we first require the following lemma. 



Lemma 3.5 For all f G C^{T$(M)), we have Ef = Ef. 

Proof. 7* = E ± f are the unique solutions to -P(7 ± ) = / with support in J^supp/) 
and 7 ± = E f are the unique solutions to -P(7 ± ) = / with support in J^supp /). Since 
trace-reversal commutes with P (Lemma 12.41) we have 



P(E ± f) = f = P(E ± f) = P(E ± f). (3.13) 
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For any k G C^iT^M)) one may show that suppA; = suppfc and hence J^suppfc) = 
J ± (suppfc). Therefore by uniqueness (from the support properties) E f = E ± f and so 
Ej=Ef. d 

Theorem 3.6 For any f G C^{S%{M)), we have Ef e ^ dD [M) if and only i/ V ■ / 6 
(n + A)C °°(7?(M)). 



Proof. By definition, Ef G S fidD {M) if and only if V • Ef = or equivalently, using 
Lemma [3.5^ V ■ Ef = 0. Taking the divergence of P(E f) = f and utilising Lemma 12.61 
we find that 

(n + A)(v-£; ± 7) = v-7 (3.14) 

and deduce that V • E ± f = E V • /, where E are the advanced and retarded Green's 
operators for (D + A) on covector fields. Hence EV ■ f = V Ef = and by [61 Thm 3.4.7], 
this holds if and only if V • / G (□ + A)C CO (T 1 °(M)). n 

We now prove an identity concerning the action of P on a pure gauge perturbation and 
then, using this, prove the relationship between E and the Lie-derivative. 

Lemma 3.7 Given a w G C°°(Tq(M)) on a cosmological vacuum background spacetime, 
then 

£(n+A)w9 = P(£ w g)- (3.15) 

Proof. Expanding out the left-hand side of (13.151) gives 

{£{n+K) w g)ab = 2V( a Dwb) + 2AV (a Wb). (3.16) 

One can show, using the Riemann tensor identity ( 11.11) . the Leibniz rule and ( 12. 2\\ , that 

V a Dw b = D(V a w b ) - AV a w b + 2R a c b d V c w d + w d V c R ach d . (3.17) 

We know from the proof of Lemma 12.61 that V a -R bcd a = on cosmological vacuum back- 
ground spacetimes and therefore we have 

V a Dw b = n(V a w b ) - AV a w b + 2R a cd W c w d . (3. 18) 

Combining this with (13.161) gives the final result. □ 

Lemma 3.8 Given a v G C^(Tq(M)) on a cosmological vacuum background spacetime, 
then 

£ kv g = E(£ v g), (3.19) 

where E are the advanced and retarded Green's operators for (□ + A) on vector fields. 



Proof. Using Lemma [321 -P^^flO = ^(n+A)^*^^ = £v 9- Thus £^± v g = E £ v g and 
the result follows by uniqueness of solutions with past /future-compact support. □ 

We are now able to prove the main result of this subsection. 
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Theorem 3.9 Any 7 G S?{M) is gauge equivalent to a Ef for some f G C^^S^M)) 
satisfying V ■ / = 0. 

Proof. Combining Corollary 12. 3[ Lemma I3T41 and Theorem 13 . 6 1 gives 7 ~ Ef with V • / = 
(□ + A)v b for some v G C™(T^(M)). Thus we also have 7 ~ Ef - £~ Ev g = E(f - £ v g), 
by Lemma 13.81 Set f '■= f — £ v g, which is smooth and compactly supported on M. 
Calculating the divergence of the trace-reversal of / gives 



V ■ / = V • / - V ■ (£ v g) = (□ + A)u b - (□ + A)^ = (3.20) 

using Lemma [2.11 □ 

The remaining lemmas of this subsection are used in section 14.41 We will require the 
notion of past/future compactness of a subset of spacetime. A subset ScMis said to be 
past /future compact if J~(p) fl S or J + {p) fl S is compact for all p G M. 

Lemma 3.10 Given f e C^(S%(M)), if 7 G ^(M) solves L ab (-f) = f ab with supp 7 
compact to the past/ future then 7 ~ — 2E f. 

Proof. By Theorem 12.21 (and the remark thereafter) there exists a«)G C°°(Tq(M)) with 
support compact to the past/future such that 7' = 7 + £ w g obeys V • 7' = 0. Using 
that L ab (£ w g) = we have L ab {-f') = L ab (j) = f ab , which simplifies, on account of the de 
Donder condition in conjunction with Theorem 12. 51 to -P(t~') = —2f . The solutions, to this 
inhomogeneous equation, with past/future compact support are 7' = —2E ± f. Lemma 1331 
entails that 7' = — 2E f. Undoing the gauge transformation gives the required result. □ 

Lemma 3.11 Given a f G C%?(S%(M)) satisfying V • / = 0, suppose that Ef = E£ v g 
for some v G Cg°(l#(M)). Then there exists h G Cg°(55(M)) such that 

f = -2L(h). (3.21) 

Proof. E(f - £ v g) = and so by [6, Thm 3.4.7] 

J=£ v g + P(h) (3.22) 

for some h G C^° (S^M)) . But for divergence- free /, an application of Lemmas 12. II and 12. 61 
gives (□ + A)(v^ + V • h) = and hence v b = —V • h. Reinserting this in ( 13. 22ft . trace- 
reversing and applying Theorem 12.51 gives the result. □ 

4 Phase Space and Quantization 

4.1 Phase Space 

We now construct the (complexified) phase space for linearized gravity on cosmological 
vacuum background spacetimes. Initially we consider the space S?(M), which by the 
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results of [35J (equation (2.21) and onwards in that reference, though our conventions differ) 
applied to the Lagrangian (I2.6p . can be endowed with a complex-bilinear pre-symplectic 
product, whose action on perturbations 7 X ,7 2 G 5^{M) is 

Ml\ 7 2 ) = / (7a6< - ll h <)dvol h , (4.1) 

where E is a spacelike Cauchy surface with future-pointing unit normal vector n, dvolh 
denotes the volume element on E associated with the induced spatial metric h and 7r is 
defined in terms of the covariant conjugate momentum II, given in (I2.10p . by 

7r afe := -n c U cab . (4.2) 

(Note that n c is pasi-pointing as a covector, owing to our signature convention.) The 
product (14.11) is independent of the choice of Cauchy surface. 

Lemma 4.1 Given 7 1 , 7 2 G ,5?(M) and two spacelike Cauchy surfaces E, E' then u^^ 1 , 7 2 ) 

Proof. Defining the current of 7 1 and 7 2 to be j c (7 1 ,7 2 ) := lab^i ab ~ ldb^-2i the P re ~ 
symplectic product of these perturbations is thus 



Ws(t\7 2 ) = / n c f(-f\~f 2 )dvol h . 



(4.3) 



Now, the divergence of the current is V c j c = / ~f 2 b L ab ( / y 1 ) — r j] lb L ab (~f 2 ) = 0, where we have 
used (I2.12p and symmetry properties of S abcd and T abcde f . Using the divergence theorem 
over the region bounded by the two Cauchy surfaces E, E' gives the desired result. □ 

Due to the preceding lemma, the E will be dropped from cue from this point on if we 
are dealing purely with solutions. 

To make the pre-symplectic product into a symplectic product, it is necessary to account 
for the degeneracies of (14.11) . that is, non-trivial solutions whose pre-symplectic product 
with all solutions is zero. The subspace of degeneracies is also known as the radical of the 
pre-symplectic form u. The next lemma shows that u is gauge invariant and therefore (as 
is well-known) pure gauge solutions are degeneracies on the space of solutions. Here, we 
work with the broader class of pure gauge solutions 5f (M) defined in (12.18D : that is, those 
spacelike-compact pure gauge solutions induced by arbitrary smooth vector fields. Recall 
that ^(M) coincides with £f (M) if M has compact Cauchy surfaces, but the two can (but 
do not always) differ in the non-compact case. 

Lemma 4.2 W(M) is contained in the radical of u. 

Proof. Suppose w G C°°(Tq(M)) and 7 G y(M), and let E be a smooth spacelike 
Cauchy surface. From Theorem IC.ll we have the identity 

We(7, £ «>9) = 2 / w b Cf{V&t^{ 1 ))dvol h (4.4) 
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and the right-hand side vanishes because C E (Datas(7)) = 0. □ 

At least in the case that M has compact Cauchy surfaces, we may prove that this 
exhausts the space of degeneracies of (j4.ip . It is natural to conjecture that the same is 
true for a large class of background spacetimes with non-compact Cauchy surfaces as well. 

Theorem 4.3 If M has compact Cauchy surfaces, the radical ofuj is precisely the subspace 
of pure gauge solutions %?(M) (which coincides with &(M) in this case). That is, given 
i G S*(M) such that w(Y,7) = for all 7 G SP{M), then 7' G &(M). 

The proof requires results of Moncrief on the ADM formulation, and is given in Appendix [Bl 
In any spacetime for which @(M ) is the radical of u, we obtain the complexified phase 
space as the quotient space 

&(M) := y(M)/W(M) (4.5) 

with weakly non-degenerate symplectic product 

"(M, [7 2 ]) = / (7>f " ll b < b )dvol h . (4.6) 

As this is independent of the choice of representative we may choose de Donder represen- 
tatives 7* (i — 1,2) for each class, for which the associated momenta are tt^ = n c @ [7J, 
where the differential operator 

@ cab Yy\ = -V c 7 afc V fe 7 ca V a 7 cfc 

has the property V c f2? mb [7] = ^(P~y) ab — h.^ ab for de Donder 7 on cosmological background 
spacetimes. 

Under complex conjugation, w([7 1 ], [7 2 ])* = a;([7 1 *], [7 2 *]). The real phase space 
&m.(M) is obtained by restricting all the above definitions to real- valued solutions and 
real-valued gauge transformations. 

We now wish to find out what the symplectic product is in terms of a de Donder 
resprentative written in terms of the solution operator, as in Theorem | 



Theorem 4.4 Given 7 e S^(M) and f G Cg°(5§(M)) satisfying V • / = 0, then 

u{[Ej\Al\) = ~\ I lTf ab dvol g , (4.7) 

z J M 

where ~f dD denotes a de Donder representative o/[7J. 

Proof. By Theorem 13.61 Ef is a de Donder solution (and its proof shows that E f also 
obey the de Donder condition). If we select a de Donder representative ~f dD of [7] then 
the left-hand side of ( I4.7J) may be written as 

u([Ej\, [7]) = / ({EJU^ ab [Y^} - Y a ^ cab [Ef]) n c dvol h (4.8) 
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using that Ef = Ef = Ef from Lemma 13.51 

As supp / is compact we may choose Cauchy surfaces E, E' such that £ C / + (£') and 
supp / C / + (E') PI /"(E). The region bounded by these two Cauchy surfaces is henceforth 
denoted by V. We will utilise the Gauss Theorem applied to the vector field 



v c = l d J?$ cab [E + f] - (E + f) ab ^ cab [Y D ] (4.9) 

on the region V to prove the desired result. Using the formula V ' c Q> cab \^j\ = ^P(~f) ab —A'y ab , 
applied to the de Donder perturbations E + f and ~y dD , the divergence of v is calculated to 

be 
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V • v = ^Tf ab + (y cl d °W ab [E + f] - (V c (E + f) ab )S> cab [Y D } (4.10) 

where we have also used the fact that P( , y dD ) = and P(E + f) = f. Using the de 
Donder condition again, the second and third terms may be seen to cancel and we have 
V • v = ^abf- By Gauss' Theorem applied to the region V, where dV = EUS' and 
with n denoting the future-pointing unit normal vector on dV, this gives 

\ [ lTf ab dvol g = - [ n c v c dvol h + / n c v c dvol h (4.11) 

^ Jv Jt, Jy,' 

The integral over £' is zero because E f and its derivative vanish on £', so we obtain 

\ J lTf ab dvol g = - j ' n c {iS&^^i] - {E + J) ab ^ ab ^}) dvol h . (4.12) 

As E~ f and its derivative vanish at E we may replace E + by — E and use (14.81) to obtain 
the final result. □ 

4.2 Observables 

Observables are functions on the (complexified) phase space, &(M). As for the scalar [TJ] 
and electromagnetic fields [TB], the observables that will form the basis for the quantum 
theory will be certain smeared fields. In our case, we wish to consider integrals of the form 
/ 1abf ab dvol g) where 7 G y(M) and / G C£°(T 2 (M)); however, this will only be gauge 
invariant, that is, independent of the choice of representative of the equivalence class of 7, 
and thus a well-defined function on ^(M), if / is restricted. 

Lemma 4.5 For f G qj°(T 2 °(M)) ; we have J M j ab f ab dvol g = for all 7 G &(M) if and 

only ifW a f\ ab ) = 0. 

Proof. For 7 = £ w g G £f(M), we have 

(V {a w b) )r b dvol g = f V a (w b fW)dvd g - f w b {V a f^)dvol g 

= - f w b (V a f {ab) )dvol g , (4.13) 
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where we moved the symmetrization to / and then used the Leibniz rule and the divergence 
theorem together with the support properties of / to obtain (j4.13p . For (14. 13[) to vanish it is 
clearly sufficient that V a f( ab ) = 0; as w may, in particular, be any element of Cq°(Tq(M)), 
necessity holds as well. □ 

We thus arrive at the final definition of observables for our theory. 

Definition 4.6 For each f G C^(T$(M)) satisfying W a f {ab ) = the observable F f : 
&(M) — > C is given by 

Ff{\ 1 }) = f lab f ab dvol g (4.14) 

(and is necessarily gauge invariant). 

The observables (I4.14p satisfy four important relations. The reader who is familiar with 
the algebraic formulation of the real scalar field might expect that we would only state (iv) 
for compactly supported /; the reason for our use of / with time-compact support will 
become clear in the next subsection. For (iv), recall that L is defined on arbitrary tensors 
by (6); in particular it vanishes on antisymmetric tensors. 

Theorem 4.7 Given any [7] G &(M), the Ff's satisfy: 

(i) Complex linearity: F a j, + gf([y]) = OiFf([y]) + f3F^([j}) for all a, (3 G C and all 

f, f G C?(T°(M)) satisfying V a f {ab) = = V a f{ ah) ; 

(11) Hermiticity: ^([7])* = F r ([7*]) for all f G C °°(7?(M)) satisfying V a f {ab) = 0; 
(Hi) Symmetry: ^([7]) = for all antisymmetric f G C^iT® (M)); 
(iv) Field equation (P22D holds: F L{f) ([j]) = for all f G C%> C {T${M)) with L(f) G 

C?(S° 2 (M)). 

Proof, (i), (ii), (iii) are obvious, and (iv) holds because L ab is formally self adjoint. □ 

We can now give a variant of Theorem 14.41 in which there is no longer any need to 
work with de Donder representatives. The theorem is proved on the assumption that 
/ is symmetric and has vanishing divergence; of course any / whose symmetric part 
is divergenceless may be decomposed into symmetric and antisymmetric parts and the 
antisymmetric part does not contribute to -P/([7J) by part (iii) of Theorem 14.71 

Theorem 4.8 Given [7] G @>(M) and f G C£°(£i?(M)) satisfying V a f ab = 0, then 

Ff(h]) = [ la b f ab dvol g = -2u([EJl [7]). (4.15) 

JM 

Proof. As V a f ab = we can use Theorem 14.41 to give 

uj([EJ], [7]) = -\j l d ab ° f ab dvol 9 . (4.16) 

As / satisfies the requirements of Lemma 14.51 we may replace 7 dD by 7 and the result 
follows. □ 
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The final two results of this section make use of the weak non-degeneracy of the sym- 
plectic product, in contrast to those above. First, we show that there are sufficiently 
many observables to distinguish points of £P(M); by this we mean that given two distinct 
equivalence classes of solutions [7] and [7'] then there exists at least one / G Co°(T 2 (M)) 
satisfying V a /(aft) = such that Ff([j}) 7^ Ff([ / y']). We are grateful to Atsushi Higuchi for 
raising this question with us. 

Theorem 4.9 Assuming weak non-degeneracy holds [thus, in particular, for any M with 
compact Cauchy surfaces] then, for any distinct [7-J, [7 2 ] G &(M), there exists a f G 
C™(S° 2 (M)) with V a /«6 = such that F f {[^ x \) ± F f (\ l2 }). 

Proof. By weak non-degeneracy there exists a [7] G £?{M) such that 

^([7],[7i])^^([7],[7 2 ]). (4.17) 

By Theorem EU [7] = [Ef] for some / G Cg°(5§(M)) satisfying V • / = 0. Using 
Theorem 14.81 together with ( I4.17P gives 

FfiM) = -M[EJ], [7i]) ^ ~M[Ej], [7 2 ]) = F f ([j 2 ]). (4.18) 

Finally, we compute the Poisson bracket of two observables in our class. Here, we 
regard &>(M) as an infinite-dimensional symplectic manifold, with the smooth structure 
determined (as a Frolicher space — see, e.g. [3H Ch. 23]) by the symplectic form u. Thus 
a curve c : R — > &(M) is defined to be smooth if t h-> uj(v, c(t)) is smooth for all v and a 
function F : &(M) — > C is defined to be smooth if F o c is smooth for every smooth curve 
c; in particular, u itself is a smooth function in each slot separately, and therefore our F/s 
are smooth. The Poisson bracket of two smooth functions F,Gg C°°(^P(M)) is given in 
terms of the exterior derivatives of F and G by 

{F,G}([<y]) = dF(IdG)\ h] , (4.19) 

where the Hamiltionian vector field IdG induced by G satisfies 

co h] (IdG\ hh v)=dG\ h] (v) (4.20) 

for v G TLi^(M) (we will show that this is uniquely defined in our context). Here ojua is 
the symplectic form at [7] G &(M). Under the identification T h] ^(M) = &(M), w [7] is 
replaced by u. 

Theorem 4.10 Assuming weak non- degeneracy holds [thus, in particular, if M has com- 
pact Cauchy surfaces], the Poisson bracket of two observables satisfying Definition ^. 0] is 

{F fl F f } = -2E{r,T)=M[Ef S \,[EF\), (4.21) 

where f s denotes the symmetric part of f [i.e., f^ b = f( a b)] and 

E(f,n := / f {ab) (EJ% d) )dvol g . (4.22) 

J M 
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Proof. We note that dFf | H ( [7'] ) = Ff ([-/]) by linearity of F f . Thus, upon using (I4.20p 
and then Theorem 14.81 we have 



u(IdF f \ h] , [7]) = *>([ 7 ]) = F r ([i\) = -M[ET\, [7]). (4.23) 

By weak non-degeneracy, this gives I dFf = — 2[Ef s ] and thus 

{Ff,F f }([y}) = -dF f \ b p\EY a \) = -2F f ([En) = -2F r {[EY s ]) 

= -2E(f s ,f s ). (4.24) 

Finally, we use Theorem 14.81 to obtain the final equality in (I4.2ip . □ 

Although superficially they do not appear the same, in fact the propagator 2E(f s , f' s ) 
is the same as the one considered by Lichnerowicz in equation (21.3) of [36]. If one expands 
out the trace-reversal then 

{F f ,F f } = -2E(f,n = -2E(f,f s ) + E(f,f), (4.25) 

where f = f a a denotes the trace of /; note the appearance of the scalar propagator E in 
the last term of this equation. 

4.3 Reformulation of the phase space 

In the next subsection we will apply Dirac quantization to the observables discussed above 
to obtain a *-algebra of observables. An alternative approach would be to directly quantize 
the complexified symplectic space £P(M) as an infinitesimal Weyl algebra. Similarly, the 
Weyl algebra provides another quantization based directly on the real symplectic space. 
(See, e.g., [19] for a presentation of both constructions emphasizing their functorial nature.) 
To clarify the relationship between these prescriptions and Dirac quantization, we now show 
that the observables form a complexified symplectic space, with their Poisson bracket as 
the symplectic product, that is symplectically isomorphic to &(M). 

Let &{M) be the set of / G C^{S\{M)) that are divergence-free, i.e., V a f ab = 0. Then, 
by Theorem 14. 7\ F : / 1— > Ff is a linear map, intertwining the complex conjugations of 
JP(M) and of complex- valued functions on <^(M) and whose range coincides with the full 
set of observables in our class. In cases where ^P(M) is weakly non-degenerate, it is clear 
from Theorem 14.81 that the kernel of the map F is precisely the subspace of / G &(M) for 
which Ef g £f (M). The next result will characterize this kernel in an attractively simple 
way. 

For this subsection only, we consider extensions of various Green's operators from 
smooth compactly supported tensor fields to tensor fields with time-compact support, i.e., 
the support lies between two Cauchy surfaces. The retarded/advanced Green's operators 
are extended by defining E f for / G Cy^Tf (M)) to be the unique solution to P( 7 ) = / 
that vanishes to the past /future of the support of /, and we define Ef = E f — E + f. 
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Many standard results have analogues for these extensions, that follow by the same argu- 
ments as the standard versions. In particular, P has trivial kernel in Cyc»(T 2 (M)), the 
kernel of the extended operator E is precisely P(C^ C (T2(M))) and any smooth solution 
to P{~y) = may be written in the form Ef for some / G C^ci^iM)). These facts may 
be summarised in the commutative diagram 

> C%> C {T$(M)) — ^-> Q? C (T 2 °(M)) — ^-> C°°(T 2 °(M)) — ^-> C°°(T 2 °(M)) 



> C -(T 2 °(M)) -^ Cg°(7j(M)) -=- ^ C (T 2 °(M)) -^ C&C3W) 

in which the upper and lower lines are exact sequences of linear maps and the vertical 
arrows are the obvious inclusion maps. The lower line is the content of [61 Thm 3.4.7]. 

Similar results hold for the extension E of the analogous vector Green's operators. 
Moreover, all the results from section |3~3| concerning the intertwining of the Green's oper- 
ators with various other operations, in particular Lemma 13. S\ carry straight through and 
apply to the extensions. 

The kernel of the map F may now be determined. In the following, Jf(M) = {L(k) : 
k G C °°(£ 2 (M))}, while JSf(M) = {L(k) : k G Cp c (S$(M))} n C °°(S°(M)). A conse- 
quence of the next result is that #(M) = &(M) if and only if J£f(M) = Jz?(M). 

Lemma_4.11 Suppose f G J^(Af). Then Ej G &(M) if and only if f G JSf(Af); simi- 
larly, Ej G &(M) if and only if f G Jgf (M). 

Proof. (=>-) Suppose Ef = £ w g G Sf(M). As V • / = 0, £ w g is a de Donder solution. 
In consequence, w satisfies (□ + A)w = and may be written as w = Ev for some 
v G Cj?c(To(M)). Therefore we have £ w g = £^ v g = E£ v g, using (the analogue of) 
Lemma I3~51 Thus Ef = E£ v g and so 

£ v g = J + P(k) (4.26) 

for some k G C^ 7 (S' 2 (M)). Taking the trace- reversal and then the divergence of this 
equation gives (n+A)tr = (D+A) V-fc, from which one can deduce, using time-compactness 
of v and V • k, that v° = V • k. Substituting this result into (I4.26P and using Theorem 12.51 
gives / = £ (y .k )i g - P{k) = 2L(k). Therefore / = L(2k) G Jz?(M) as required. 

(<=) Conversely, given k G C%P C (S%(M)) satisfying L{2k) = J, let w = Eiy ■ fc) s G 
C°°(Tq(M)), which gives £ w g = £e(v-W^ = -^Av ■&)»#' us i n § (the analogue of) Lemma [3781 
in the final equality. Now by Theorem 12.51 this becomes £ w g = E(2L(k) + P(k)) = Ef, 
because EP{k) =0. As / is compactly supported by assumption, we deduce that 
£ w g G C^ c (s{(M)) and hence Ej G #(M). 

The second statement has an exactly analogous proof, replacing C™ c by C^° and C°° 
by CgQ, and hatted spaces by their unhatted counterparts throughout. □ 
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It follows that the vector space of observables is isomorphic to J£~(M)/J2f (M) if weak 
non- degeneracy holds. 

Proposition 4.12 The map f+J&(M) !->■ [Ef] is a linear isomorphism of &{M)/J£(M) 

and 3P{M\ Accordingly, if weak non- degeneracy holds [thus, in particular, for any M with 
compact Cauchy surfaces], there is an isomorphism of the space of observables and &(M) 
such that Ff h-> [Ef] for all f G &(M), and which is symplectic if the observables are 
equipped with symplectic product a(Ff 1 ,Ff 2 ) = j{Ff 1 ,Ff 2 }. 

Proof. Lemma [A. Ill shows that the given map is well-defined and injective, while surjectiv- 
ity is evident by Theorem 13. 91 The symplectic product is seen to be preserved by combining 
Theorems KS\ and KM to give 4w([-E/ 1 ], [Ej 2 ]) = {F fl ,F f2 } for any f x , f 2 G &{M). u 

4.4 Quantization 

To quantize the theory we follow Dirac's prescription, which requires that we seek operators 
Ff, labelled by test tensors / with divergence- free symmetric part and whose commutator 
is given by [Ff, Ff] = i{Ff, Ff}l. We also expect these operators to respect the quantum 
analogues of the relations for the classical observables from Theorem 14.71 As we will 
regard these objects as smeared quantum fields, we use the notation [7J(/) rather than 
Ff, cautioning that the [7] in such expressions is not to be confused with an equivalence 
class of classical solutions. Combining the above requirements gives the algebraic relations: 
(i) Complex-linearity: [y]{ocf 1 + f3f 2 ) = «[7](/i) + P[l](f2) f° r all a,/3 G C and all 
fi G C™(T°(M)) such that V a (f) (ab) = 0; 

(ii) Hermiticity: [t](/)* = [7K/* ) for all / G C^°{T^(M)) such that V a f {ab) = 0; 

(iii) Symmetry: [y](f) = for all antisymmetric / G C^°(T2 (M)); 

(iv) Field equation: [~f](L(f)) = for all / G C^ C (T^(M)) such that L{f) G C °°(^(M)). 
(v) Commutation relation: [[ 7 ](/i), [7K/2)] = -^E{f{, f s 2 )t for all /, G C^{T^{M)) 
such that V a {f) {ab) = 0. 
Note that (a) in the statement of the field equation, the linearized Einstein tensor L(f) 
is symmetric (by definition of L) and divergence free (by the Bianchi identities) for any 
/ G Ctq(T 2 (M)); (b) the final relation implies that the commutator vanishes if the f i are 
spacelike separated, as required by the Bose statistics of the spin-two field. As discussed 
in section I4.3[ the use of time-compact test tensors in (iv) is necessary to allow for certain 
effects of nontrivial topology. An analogous modification to the axioms for the electromag- 
netic field is required if the topological restrictions imposed in [TBI E] are relaxed. Also 
note that (ii) corresponds to item (ii) of Theorem 14.71 for real- valued 7. 

The algebra of observables A(M,g), where (M,g) is the background spacetime, will 
consist of finite linear combinations of finite products of [7K/), [7](/)* and a unit 1 
obeying the above relations. More formally, it may be constructed as follows: we first form 
the free unital *-algebra generated by symbols [7K/) together with a unit 1. However, 
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this algebra is too large as it does not take into account the relations. To impose them we 
quotient by the two-sided *-ideal they generate 

The algebra A(M, g) may also be constructed from a different perspective, as the 
application of a quantization functor to the classical complexified phase space £P(M) with 
symplectic product to (see [191 Sec. 5] for the definition and properties of this functor). 
This is particularly convenient in the discussion of covariance below. The equivalence of 
these approaches follows directly from Proposition 14. 121 

We now state and prove a time-slice property. This shows that the algebra is generated 
by smeared fields with smearings supported in any slice around a Cauchy surface. 

Theorem 4.13 Given any connected causally convex open neighbourhood JV of a spacelike 
Cauchy surface £ and a f e C^(S$(M)) with V-/ = then there exists afe C^C^MO) 
with V • / = and ahe C%°(S$(M)) such that 

f=~f-2L(h). (4.27) 

This entails that [j](f) = [j](f) in A(M,g). 

Proof. This follows the method used in the electromagnetic case from [TT| Prop. A. 3(b)]. 
By assumption, jY is globally hyperbolic, and we may choose Cauchy surfaces S 1 * 1 for M 
with Y^ 1 C / ± (S) n JY . Take two scalar functions x ± £ C°°(M) satisfying x + = 1 in 
J+(£+), x + = in J-(E-) and X + + X~ = 1- Define 

f:=2L( X + Ej), (4.28) 

which satisfies V ■ / = by the linearized Bianchi identity and is compactly supported 
within jY (x + Ef evidently vanishes to the past of JV and coincides with a de Donder 
solution to the linearized equations to the future of JV by hypothesis on /). Note that 

( 14728]) implies that 2L(x~ Ef) = -f. By Lemma EHH] we have -E + f ~ x + Ej, E f ~ 
X~Ef and hence 

Ef = Ej+£ w g (4.29) 

for some £ w g 6 £f (M). As V ■ / = V ■ / = 0, both Ef and Ef are de Donder solutions 
and so w solves (□ + A)w = (see the remarks following Theorem I2.2f) : hence by [HI 
Thm 3.4.7], w = Ev for some v e C^(T^(M)). Substituting this result into ( 14T291 and 

using Lemma I3T51 gives E(f — f — £ v g) = 0. Using Lemma 13.111 gives the desired resultn 

Finally, we note that the theory could also be quantized by means of the Weyl algebra. 
This results in a unital C*-algebra W(M,g) generated by elements WQ7]) ([7] G ^ R (M)) 
obeying the relations 

W{0) = 1, W([ 7 ])* = W{-[i\), W([i\ + [V]) = e^(W'M)/ 2 W([7])W([7l); 

and with the unique norm making the resulting system a C*-algebra - for an explicit 
construction of the algebra see, for example, [37] or [B] (example 4.2.2 of that reference) in 
the case of non-degenerate u, or [9j in the degenerate case. 



d We use the same notation for [-f](f) and the corresponding equivalence class in the quotient. 
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4.5 Covariance 

We briefly discuss the extent to which the theory we have constructed may be formulated 
as a locally covariant theory in the functorial sense introduced by Brunetti, Fredenhagen 
and Verch (BFV) [TT]; however, for brevity, we will not emphasize the categorical struc- 
tures here (they are easily reinserted). Let (M,g M ) and (N,g N ) be globally hyperbolic 
spacetimes solving the vacuum Einstein equation with cosmological constant A, and as- 
sume that these spacetimes are endowed with time-orientations and orientations, which we 
leave implicit in what follows. Consider any smooth embedding ip : M — >■ N that is an 
isometry, preserves the (time-)orientation and has a causally convex image. We restrict 
to spacetimes for which &(M) and &(N) are both weakly non-degenerate, though we do 
not assume that M and N have compact Cauchy surfacesjj 

Owing to the isomorphism in Proposition I4.12[ the push-forward ip* of compactly sup- 
ported tensor fields induces a linear map &{ip) : &{M) — > £P(N) so that £P(ip)[E M f] = 
[E N ip*f], provided that ip*Jf(M) C Jt?(N). This raises the following question: given any 
k G C^ 7 (S , 2(M)) so that Ljw(fc) is compactly supported, does there exist k' e Crc^l^O) 
so that Ljv(fe') = ip*L M (k)l Clearly, if k is compactly supported, k' = ip*k has the re- 
quired property, so ip*J£(M) C Sf(N) C Jzf(iV). However, this argument cannot be used 
if A; is a general element of Cj^S^M)), and it seems that there may be a genuine ob- 
struction in some (though not all) cases where Jzf (M) is a proper subspace of Jzf (M). This 
is analogous to the distinction between de Rahm cohomology with or without compact 
support; an equivalent question here is whether or not ip*@(N) exhausts £f (M). 

If, indeed, ?/>*Jzf (M) C Jzf (iV) then the map &(ip) is easily seen to be symplectic, as 

u N (^(iP)[E M J], @>{il))[E M J]) = u N ([E N ipJ], [E^J]) 

= - l -E N {^f^J) = - l -E M {f,J) 

= ou M ([E M J},[E M f}) (4.30) 

for arbitrary /, / e JzF(M). It is therefore also injective because &(M) was assumed to 
be weakly non-degenerate. It is also clear that &(ip) commutes with complex conjuga- 
tion. Consequently, the general properties of the quantization functor used to construct 
A(M,qm) from 3?{M) (see [19]) entail the existence of an injective, unit-preserving *- 
homomorphism A(ip) : A(M, g M ) — >■ A(N,g N ), uniquely determined by the property 
«4-(V')[7]m(/) = [t]v(V ; */)- If we consider an additional embedding ip : M' — > M of the 
above type, it is evident that £/(ip) o sf(ip) = s$(ip o tp); it is also clear that srf maps the 
identity mapping of M to the identity mapping of &/(M,g M ). Moreover, if the image of 
ip contains a Cauchy surface of N, then the map £/(ip) is surjective by Theorem 14.131 and 
hence an isomorphism. 

This discussion may be summarised in the language of [UJ by saying that stf defines a 
covariant functor from a certain category of globally hyperbolic spacetimes to the category 



c Note, however, that if M has compact Cauchy surface, the existence of the embedding ip entails that N 
also has compact Cauchy surfaces, which are oriented-diffcomorphic to those of M; see [18j Prop. 2.3(a)]. 
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of unital *-algebras with unit-preserving injective *-homomorphisms as morphisms, and 
that this theory has the time-slice property. Similarly, the quantization of ^r(M) in 
terms of Weyl algebras defines a functor to the category of unital C*-algebras. However, 
the category of spacetimes is a subcategory of that usually studied in the BFV formalism: 
the spacetimes themselves are restricted to those cosmological vacuum solutions on which 
the phase space of linearized gravity is weakly non-degenerate, while the morphisms must 
be restricted to permit only those embeddings ip : M —¥ N for which ^^{M) C Jzf(iV). 
A more geometrical characterization of this class would be desirable. 

The restriction to cosmological vacuum solutions has an important consequence: it is 
not possible to formulate the relative Cauchy evolution that plays an important part in 
the BFV formalism - see [11] and [T8| [T9] - but which requires the freedom to consider 
arbitrary (sufficiently small) compactly supported metric perturbations. This is related to 
the well-known absence of a local stress-energy tensor for the gravitational field. Finally, 
the detailed study of circumstances under which &(M) is weakly non-degenerate, and of 
those morphisms ip '■ M — >■ N for which &*{$) is well defined, are evidently interesting 
problems with a cohomological flavour and deserve further investigation. 
Acknowledgements: The authors thank Atsushi Higuchi for many extremely useful dis- 
cussions on this subject and Vincent Moncrief for invaluable assistance (via Atsushi) in 
understanding non-degeneracy of the symplectic product. We also thank Ian Mcintosh for 
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A Transverse-traceless gauge 

A.l Conventions 

As the proof of Theorem 12.71 makes use of differential forms we briefly summarize the 

conventions employed. We work with the smooth spacetime (M,g) satisfying all of the 

topological criteria from section [2j Our conventions follow [lj and are consistent with 

those of [T71 H3]: for the space of p-forms on M we use the notation Q P (M), the space of 

compactly supported p- forms on M is denoted fig(M) and the space of p- forms that are 

spacelike-compact is denoted by Q P SC (M). 

The wedge-product of ex E W(M) and f3 E Q q (M) is ex A (3 E fi p+9 (M), which is given 

by 

(p -\- q\\ 

[a A p)a 1 ...a p+q — -j-j <X[a 1 ...a p Pa p+ i...a p+q ]- (A.l) 

The Hodge star operator is the map * : Q P (M) — > Q n ~ p (M) uniquely defined [TJ 
Prop. 6.2.12] by the condition ex A */3 = (ex, (3) g dvol g for all ex,f3E Q P (M). Its square is, 
by [U Prop. 6.2.13], equal to 

(*) 2 = (_i)p(»-p)+» > (A.2) 
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where n is the dimension of the manifold and s is the index of the metric g. On a four- 
dimensional Lorentzian manifold (in — h++ signature) this becomes (*) 2 = (— l) p+1 , while 
on a spacelike Cauchy surface thereof it reduces to (*) 2 = 1. 

There is a standard pairing, see (TJ p. 538], between p-forms on a manifold M: given 
w G VP(M) and / G fig(M) one defines 



(w,f)u--= / wA*f. (A.3) 

J M 

If M is compact then the restriction that / be compactly supported can be dropped. 
The exterior derivative d : Q P (M) — > Q P+1 (M) is given by 

(daOai-op+i = (p + l)V [ai a a2 ... ap+l] , (A.4) 

where a G fiP(Af). The codifferential <5 : QP(M) -> Qp-\M) is, by [1, Dfn 6.5.21], given 
by 

( j = (_l)»(p-i)+»+i*d* (A.5) 

and always annihilates 0-forms. On the spacetime (M,g), 5 = *d*, while on spacelike 
Cauchy surfaces, 5 = (— l) p * d*. 

By using Stokes' Theorem one can show that 

(da, 13) = (a, 8(3} (A.6) 

for a G f2 p_1 (M) and /3 G fi p (M), provided at least one of them is compactly supported. 

A. 2 The transverse-traceless gauge 

By using differential forms we can make use of the methodology of [US Prop. 2.6], which 
deals with obtaining the Lorenz gauge in electromagnetism, where there is an identical 
hyperbolic equation and a constraint similar to (12.411) . 

As in [4~3l Sec. 2.4], define % : £ — > M to be the embedding of the Cauchy surface S in 
the spacetime M and define the following forms on a Cauchy surface E: 

(A.7) 

(A.8) 

(A.9) 

(A.10) 

which together constitute the Cauchy data on S for w, meaning they correspond to iu|s 
and V n w\z. We sometimes use the notation p( ), P(d), P{5) and p( n ) from J43J Sec. 2.4] for 
the linear maps (1A.7I) . (1A.8I) . (1A.9J) and (1A.10I) respectively, applied to one- forms and zero- 



W (0) 


= i*w G nl(E) 


w (d) 


= -*i* *dw G f^(£ 


W {5 ) 


= i* Sw G fig(E) 


U)( n ) 


= -*r *w eng(E), 



forms. The forms on a Cauchy surface corresponding to the Cauchy data for a zero-form 

7 are 

7(o):=P(o)7e^o(S) (A.ll) 

7(d):=P(d)7efi°(S). (A.12) 
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(Note: 7(5) and 7( n ) are automatically zero.) Writing (j2.4ip in forms notation gives 

-(Sd + dS)w + 2Aw = 0, (A.13) 

Sw = -7 (A. 14) 

respectively. In (IA.14I) we used 5w = —V a w a and in (1A.13I) the 'extra' Aw comes from 
non-commutativity of covariant derivatives and using (12. 2\\ . 

We will require scalar and one-form Green's identities, see [281 Appx A], [131 Sec. 2 
which link solutions of hyperbolic equations to their initial data on a spacelike Cauchy 
surface E. For w G Qg C (M) solving (IA.13I) we have 

(to, f)u = {W(0),P(d)Ef)v + (W(S),P(n)Ef)x - (W(d), P(0)-E/>23 - (W(n),P(8)Ef)s, (A.15) 

where / G Qq(M) and E is the advanced- minus- retarded solution operator for the differ- 
ential operator — (Sd + dS) + 2A acting on 1-forms (see section [331 for further details on 
Green's operators). The scalar case is (1A.23I) . 

We now show what constraints the Cauchy data need to satisfy in order that the solution 
to (1A.13P also satisfies (1A.14D : we adapt [121 Prop. 2.6] to achieve this. 

Theorem A.l Suppose w G Qg C (M) solves (— (Sd+d5)+2A)w = and^f is a de Bonder 
solution, then 5w = I7 if and only ifw^g) — §7(o) an d Sw^d) + 2Au>( n .) = |7(d). 

Proof. (=>) The pull-back of the constraint 5w = ^7 to the Cauchy surface gives w^g) — 
§7(o) j while applying — * i* * d to the constraint gives 



Using (IA.13I) . this yields 



1 1 

*i* * d5w = p {d) 5w = -p(d)l = 777(d)- (A.16) 



* i* * 5dw — 2A * i* * w = -7(d)? (A- 17) 



which may be rewritten as 5w(d) + 2Aw( n ) = 57(d) by using the definitions of the various 
quantities and the fact that d commutes with i*. 

(<=) To prove that such a w will satisfy the constraint 5w = ^7 globally on (M, g) we 
begin by taking an arbitrary / G Qq(M) and computing 

(SwJ)m = (w,df) M = (w(p),P(d)Edf)v + (w { s),P(n)Edf)x 

- (w( d ), P(o)E<±fh ~ (w(n),P(6)Edfh, (A.18) 



where we use (1A.6|) and (1A.15J) 



f The formulae in [28, 43] are modified by changing to our — h++ signature convention. Note that there 
are sign errors in Eq. (2.21) in [33], and that [35J uses the retarded-minus-advanced propagator, whereas 
we use the advanced-minus-retarded. 
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From [43, Prop. 2.1] we know that -Ed = dE, where E is the advanced-minus-retarded 
solution operator for the differential operator — (<5d + dS) + 2A = — <5d + 2A acting on 
0-forms. Using this and p(d)d = 0, ( 1A.18I) reduces to 

(5w,f) M = (w {5 ),p( n )dEf)x - (w (d) ,p {0) dEf)x - (w (n) , p^dEf)^. (A.19) 

For the second term on the right-hand side we can use that the pullback and the exterior 
derivative commute, and then ( 1A.6I) to obtain 

(SwJ)m = (w( 5 ),P(n)dEf)x - (6w( d ),P(o)Ef)x - (w( n) ,p^dEf} x (A.20) 

1 1 

= 2^(o)'P(«) d ^/)s - ^(l{d),P{0)Ef)z 

+ 2A(iw (n) , P(p)Ef) - (w {n) , P{S) dEf)x, (A.21) 

where in (1A.21I) we substituted the restrictions on the Cauchy data. The last two terms 
cancel because p^dEf = p^SdEf = 2Ap^Ef. Therefore 

(SwJ)m = ^{l(o),p(d)Ef}^ - -(j(d),P(o)Ef) x , (A.22) 

where we have used that, in this case, p(d) = P(n)d. 

The trace of a de Donder solution 7 satisfies the scalar wave equation (I2.40p . which 
in forms notation is — Sdj + 2A7 = 0. Therefore we can use the scalar version of ( 1A.15I) 
[recall that 7 ((5) = 7 (n) = 0], i.e., 

(7, Dm = (7(0), P(,d)EDx ~ (7(d), P(o)E/) s (A.23) 

for any / G fig(M). Comparing this with (1A.22I) gives (5w,Dm = \{i,Dm for all / G 
Qq(M), and hence 5w = ^7. n 

We may now give the proof of the splitting y(M) = y TT (M) + <g(M) for A/0. 

Proof of Theorem [Ml We know from Corollary O that S*(M) = y dD (M) + &(M), 
therefore if we can decompose the space of de Donder solutions as y dD (M) = y TT (M) + 
Sf (M) D y dD (M) then we can achieve (T2T421 . Given a perturbation 7 G y dD (M) on a 
cosmological vacuum spacetime (M,g) with A ^ 0, the constraints of Theorem IA.1I are 
satisfied by iu (0 ) = 4^7(0), iU( d ) = 0, W( n ) = ix7(d) and w^s) = §7(o) as Cauchy data. 
The solution with this data is w = Trdj, which corresponds to the choice in equation (9) 
of [32] for de Sitter spacetime. Therefore appropriate Cauchy data exists and one may 
gauge transform from the de Donder gauge to the transverse-traceless gauge. □ 

In the case A = 0, the second constraint of Theorem IA.1I reduces to Sw^ = |7( d ), 
which becomes a cohomo logical problem. The scalar 7( d ) is co-closed, Sj^) — 0, but is it 
co-exact? Equivalent ly, we must solve 

d(*ttf(d)) = -- *7 (d ), (A.24) 
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in which *7(d) is a 3-form on E and necessarily closed. 

There are two cases to consider: depending on whether or not w has compact support 
on Cauchy surfaces. If as we assume, w has compact support on Cauchy surfaces then 
from [U Thm 7.5.19(i)], *7(d) is exact if and only if 



*7(d) = 0. (A.25) 

s 

If, on the other hand, E is non-compact and w is allowed to have non-compact support on 
E, [TJ Thm 7.5.19(iii)] gives that if 3 (E) = and so *7(d) is exact and the TT gauge may 
be attained as is standard, e.g., in Minkowski spacetime. 

B Non-degeneracy 

B.l Background on the ADM formalism 

In order to prove Theorem 14.31 we need to appeal to the results of [391 [22], which use the 
ADM formalism (for original references, see [3]). This formalism puts general relativity into 
the form of a dynamical system, where given a three-dimensional smooth manifold E and 
Cauchy data (h,va) e C°°(S%(Y,)) x C°°(S 2 (E)) we can obtain a solution ((-e,e) x E,g) 
to Einstein's equation. Here S'q(E) is the space of smooth second rank contravariant tensor 
densities on E, h is the spatial metric on E and zu ab = yh (k ab — ^h ab k) is the conjugate 
momentum, where k is the desired extrinsic curvature^ of E in the solution spacetime and 
h is the determinant of the metric h. In fact, to obtain a solution, one also needs to specify 
a lapse function and a shift vector field on E, which can both be time-dependent and are 
freely specifiable and non-dynamical. Together they make up the components of a vector 
field whose integral curves provide a flow of time in spacetime. The spacetime metric g is 
constructed from the lapse and shift as well as the evolved spatial metric h obtained from 
solving the ADM equations (IB.17P below. 

The initial data (h, vj) are not freely specifiable as they need to satisfy constraints 
given by the map $ : C°°(S%(Y,)) x C°°(S%(E)) -► C°°(E) x C°°(T$(Yi)), where 

$(h, vj) = (Je(h 7 vj), 8{h, w)) (B.l) 

and the Hamiltonian and momentum constraints are 

ab fr Jl 

JT{h, vj) = -R^\h) + ^p - |- + 2A (B.2) 

and 

S%h,V7) = D b f^\ (B.3) 



s Our convention differs from 1221; however, the overall definitions of vj coincide. 
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respectively. Here R^(h) is the Ricci scalar for the metric h and D a is the covariant 
derivative associated with h. Vanishing of (IB.lj) is a necessary condition for a spacetime 
to be a solution to the vacuum Einstein equation with cosmological constant. 

For linearized gravity in the ADM formalism, one considers a one-parameter family of 
Cauchy data (h(X), "^(A)) and takes the derivative at A = 0. Thus the Cauchy data for the 

linearized ADM equations are (7^ 3 \p) = ( dx , q X ) and the linearized constraints 

are the components of the derivative of the constraint map ( IB. II) at (h, vj). One should note 
that to solve the linearized ADM equations one needs to specify a linearized lapse function 
and linearized shift vector field; as in the full non-linear case they are non-dynamical and 
freely specifiable. 

From now on we assume that the background is a solution to the vacuum Einstein 
equation with cosmological constant and so $(/i, vj) = 0. The linearized constraints are 
the derivative of the constraint map (IB.ip . 



D${h,vj) : C°°(5 2 °(S)) x C°°(S 2 (S)) -)■ C°°(£) x r^p)) 
evaluated at (h,vj), where 



(B.4) 



D${H,vj) 



7 



(3) 



P 



(3) 



V 



D8{h,vj) 



T 
P 



J 



(B.5) 



whose actions are given by 

,(3>\ 1 



D3e{h,vv) 



T 
P 



h 



„ab^ L 2 \ „,(3) 

^7 ^ab - -K> ) 7 



+2 I zu ac w 



cb 



:^ b 1 1% 



: 2 | W^P^ ~ -UJP 



~ ( D a D b lS 



D a D al ^ - R {3)ab ^A (B.6) 



and 



D6{h,vj) 



7 



(3) 



P 



1 



2D b p ao + ti7 6c (D c7 ( 6 3)a + J D fe 7 



(3)o 
c 



^ a 7i 



(B.7) 



where 7 < - 3 - ) = h ab x ah , zu = h a bzu ab and p = h a bp ab . To get these into a form analogous to 



that in [391 122] we evaluate f ]B.6j) on the constraint surface to give 



D3tf(h,vj) 



7 



(3) 



P 



1 

h 



W ZU ab 



-ZU 2 ) 7^ 



2 w ab p' 



(lb 



1 



07£> 



+2 w ac w r 



\^ b )l^ 



D a D b 1 



ah 



2 



R (3)ab _ _ h ab R (3) + Ah ab 



(3) 
Tab 



. (B.8) 
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Note that the difference between this and the A = case considered in equation (2.8) of 
[39] is the cosmological constant term. 

We also require the following inner products, defined in [39] (equations (2.4) and (2.6) 
in that reference). The first acts on the domain of D<$>(h,vj), i.e., the vector space 
C°°(5 2 °(S)) x C°°(5 2 (S)) and is 

((7 (3) ;p) . tf$ )P)) ;= jf (Vh 1 ^Sh ac h bd + -Lh ac h bd p ab pA , (B.9) 

where ^y( 3 \^ 3 > e C°°(S , 2(S)) and p,p E C°°(S'q(S)). The second acts on the codomain 
C°°(E) x C°°(T£(E)) of D$(h,va) as follows 



(((f } V);(f,V)}}:= / (/ • / + h ab V a V b )dvol h 



(B.10) 



where /, / e C°°(E) and V, V G C^T^E)). 

The adjoint of the differential operator DQ(h, vj) is calculated, using the inner products 
(1B.9J) and (IB.lOj) and integration by parts, to be 






where DJ?(h,vp)*{f) = (a, (3) with 



(B.ll) 



"ab 



-- I G7 Ca W cd - -W Z \ h ab f +2 I tZ7 ac CJ7 C fe - -G7 afc G7 ) / 



D a D b f - h ab D c DJ - (r® - h ab R^ + Ah ab \ f 



and 



Also 



»'"'' = 2/ ( w a6 - -ro/i a6 



-Vh(D a V b + D h V a ) 



(B.12) 
(B.13) 
(B.14) 



We now define a unitary U : C°°(S2(£)) x C°°(Sg(E)) -)• C°°(5 2 (E)) x C°°(S 2 (E)) by 



t/(7 (3) ,P):= 



v 7 /^ 



=p bb ,v^( 



so that £/ o D§(h, vj)* corresponds to 7(fo, v&) 



7 



1 



(3)\M 



(B.15) 



/}$(fa, "CcO"!" from [39], where 



DQ(h, vjy is the 'new form of the adjoint' defined in equation (4.2) of that reference. The 
inverse map U' 1 : C°°(S§(E)) x C°°(S 2 (E)) -»- C°°(S 2 (E)) x C°°(S 2 (E)) is given by 



f/- 1 (7 (3) ,p)=(-^P bb ,-v / ^(7 (3) ) Stt ). 



(B.16) 
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The ADM evolution equations may be written as follows 

U«® =u - loD * {h '» ) i-»)' (Ri7) 

where N is the lapse function and N is the shift vector field associated with the slicing. 
The ADM symplectic product on the background (h,Tu) is, see [221 P- 333], given by 

<^((7 {3) ,P); (7 (3) ,P)) = / (l%P ab - l$?P ab )d 3 x. (B.18) 



Observe that 



<^((7 (3) ,p);(7 (3) ,p)) = ((7 (3) ,p);^ 1 (7 (3) ,p))). (B.19) 



B.2 Analogues of Moncrief ' s splitting theorems 

In [39] it is shown that the space of initial data can be decomposed into orthogonal sub- 
spaces, using the inner product (1B.9[) . Here, we generalize these decompositions to the 
case of nonzero cosmo logical constant. The first splitting is 

C°°(S 2 °(S)) x C°°(5 2 (S)) = ker £>$(/*, w)0 range D®(h,<&)*, (B.20) 

where ker D$>(h, vj) is the subspace of data satisfying the linearized constraints and 
range DQ(h, vj)* is the unphysical data. This splitting was done for the case of A = in 
[39\ Sec. 3] using ellipticity of the operator DQ{h,vj) o D<&*(h,z&), which is proven by 
showing that D<&*(h,vj) has injective principal symbol and applying [7, Thm 4.4], which 
is valid on compact Riemannian manifolds. Since our modifications to the linearized con- 
straint map and its adjoint only add a Ah term, the principal symbol will be unaffected 
and so the operator is still elliptic. Hence the remainder of the Moncrief argument of [3"9l 
Sec. 3] remains valid and the first splitting holds for general A. 

The second splitting decomposes the subspace ker D$>(h, vj) into a pure gauge sub- 
space, meaning data for pure gauge solutions, and a physical subspace. In [UH Sec. IV] it 
is shown that data for a pure gauge solution to the linearized equations corresponding to 
£ w g, on a vacuum spacetime with A = 0, is given by 

(7 (3) ,p) g au g c = UoD$(h,^y 0T) , (B.21) 

where C = n a w a and X a = q a b w b are respectively the normal (with respect to the future 
pointing normal vector n) and tangential projections, relative to £ (using the associated 
projection tensor q a b ), of the gauge vector field. The above result was initially proved via 
a lengthy calculation, and later by more geometrical methods, see [221 Thm 4.7]. The 
result (IB.21I) also holds on vacuum spacetimes with non-vanishing cosmological constant 
by following the same argument used in the proof of [22^ Thm 4.7] but instead using the 
vacuum ADM equations with cosmological constant (IB. 171) . 
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Before performing the final split, one needs to check that the pure gauge subspace lies 
in the constraint subspace. Again, one could check this by lengthy calculation, as was done 
in J39j Thm 4.1] for the A = case; instead, we appeal to the geometrical method of [22} 
Prop. 3.2] whose result is unaffected by the inclusion of a cosmological constant. 

With the two preceding results and, as argued earlier, ellipticity of DQ(h, vj)oD§*{h, vj) 
unaffected by addition of a cosmological constant, the subspace ker D${h, xxj) can be de- 
composed, see the argument in between Theorem 4.1 and Theorem 4.2 of [39], into 

ker D$(h, vj) = range(£7 o D$(h, vj)*) © ker((C7 o D$(h, vj)*)* n ker D$(h, vj), (B.22) 

where the first space is pure gauge and the second space is the physical space. 
Therefore the final split of the initial data is 

C°°(5 2 (S)) x C°°(5 2 (£)) = range D<S>{h,vj)* ®range{U o D<S>{h,vj)*) 

© ker((C7 o D$(h, vj)*)* n ker D$(h, vj), (B.23) 

which is the same result as the A = case from [39l Thm 4.2]. This decomposition allows 
one to prove that on the space of initial data obeying the constraints, ker DQ(h, vj), 
the only degeneracies of the ADM symplectic product are pure gauge. We now give the 
analogue of [21, Prop. 4.38]. 

Theorem B.l The ADM symplectic orthogonal complement to the subspace ker D§(h, vj) 
is the pure gauge space range(£7 o D§(h, vj)*). 

Proof. Let (7 (3) ,p) G ker D$(h, vj) satisfy w^( (7 (3) , p); (7 (3) ,p)) = for all (7 (3) ,p) G 
kerD$(h,vj). Then by flBTTOl) . 

(( 7 (3) ,P);^ 1 (7 (3) ,P))) = (B.24) 



and so U 1 ('J ,p) is orthogonal to ker D&(h,vj). By the first Moncrief split ( IB. 201) this 
means that £/~ 1 (7 ( - 3 - ) ,p) G range D§(h, vj)*. Hence 

(7 (3) ,p) Grange([/oL)$(^,-cc7)*) (B.25) 

and is therefore pure gauge. □ 

B.3 Proof of Theorem [4^1 

The main issue is to translate Theorem IB. II into the setting studied in the main body of 
the paper. Begin by taking an arbitrary smooth spacelike Cauchy surface S and denote 
by jV a normal neighbourhood of E. (For details about normal neighbourhoods, see the 
second paragraph of section [2731 ) Assume that a solution j' G S?{M) is a degeneracy of 
the symplectic form u, i.e., a; (7', 7) = for all 7 G 5^[M). Without loss of generality, 7' 
may be chosen synchronous near S; it will be enough to restrict attention to synchronous 7 
as well. (Theorem 12.81 means that we can gauge transform any solution to the synchronous 
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gauge near £ and since, by Lemma I4.2[ pure gauge is a degeneracy then 7' will still be a 
degeneracy of u>.) 

We now restrict our attention to the normal neighbourhood JY ', on which we can 
introduce Gaussian normal coordinates. In such coordinates the spacetime metric takes 
the form g = —dt <S> dt + h^dx 1 (g) dx^ and the synchronous condition is precisely 7o M = 0. 
The solutions 7', 7 correspond to solutions to the linearized ADM equations about the 
background (jV \g\j/) in the slicing given by the Gaussian normal coordinates: thus we 
have unit lapse, vanishing shift (and vanishing linearizations thereof). The corresponding 
ADM Cauchy data are (7 /{3) ,p')> (7 (3) >P) e £°°(S 2 (£)) x C°°(S$(Ti)) respectively, where 
in these coordinates 

lf=Hi\x (B.26) 

pfc = ^fh— (h im h jn - ti j h mn ) d h mn 

~ ^ (tS #" + ^ m 7S - 7g)^ n - ^7© ) do/w 
+ ^ (^m^» _ fc « r ) f(V n7 ) mn | s + ^ (ft^S + Q^tS)) ■ (B.27) 
Using, for convenience, Gaussian normal coordinates one may sIiowlj that 

W (Y,7) = ^4((V (3) ,pO; (7 (3) P))- (B.28) 

By Theorem IB . 1 1 and degeneracy of 7', (7^ 3 - ) , p') is data for a pure gauge solution; therefore 
on the region jV , 7' = £ w g for some ii> G C oo (T 1 ( I y^ / ')). Now perform a global gauge 
transformation on 7' using a vector field i> G C°°(Tq (M)) satisfying v = — iu on an open 
neighbourhood of S within jY '. The result will still be both a solution and a degeneracy in 
<y(M) but has Datas(7' — £v9) = (0, 0) and therefore by Theorem 13. 3[ 7' = £ u g for some 
u G C°°(Tq(M)). Note that due to the compactness of E, all three vector fields w, v and 
u will be spacelike-compact and hence so will their associated pure gauge perturbation. 
Therefore 7' G S?(M) = #(M). 

C Identity connecting the symplectic product and lin- 
earized constraints 

Theorem C.l For any 7 G 3T{M), w G C°°(Tq(M)) and smooth spacelike Cauchy sur- 
face £ u>z£/i nnzi future-pointing normal vector n, we have 



u*(l,£ w g) = 2 / w a L ab (j)n b dvol h = 2 / w 6 C 6 s (Data E (7))^o/ h . 



(C.l) 



h We caution that the relationship between ir and p is not straightforward, although they do coincide 
on constant time hypersurfaces in Minkowski space. 
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Proof. Choose a vector field v G C°°(Tq(M)) such that v = w in a neighbourhood of E 
and v vanishes to the far past of S. (This trick originates, as far as we know, from |27j . 
see the paragraph preceding equation (79) in that reference). By ( 13.21) we have 

w a C*{T)&%^{l))dvol h = / w a L ab (j)n b dvol h . (C.2) 



The proof now uses two identities to successively re-express the right-hand side of (IC.2J) . 
Lemma C.2 VFrf/i 7, t/; and u as above, 

[ w a L ab (-y)n h dvol h = - I V (a ^)L a6 (7)aW 9 , (C.3) 

where M~ = /~(£) is the region to the past of the Cauchy surface E. 

Proof. Using the properties of v and the Gauss Theorem on the region M~ we have 

/ w^^^dvoln = - [ V b (v a L ab (<y))dvol g} (C.4) 

where we have used that v vanishes to the far past. The right-hand side can be rearranged 
using the Leibniz rule, symmetry of L ab and V a L ab = to give the result. □ 



Now utilise the pre-symplectic product (14. lj) to re-express the right-hand side of ( 1C.3J) . 



Lemma C.3 With 7, w and v as above, 

oo s (j,£ w g) = - [ 2V {a v b) L ab ( 1 )dvol g . (C.5) 

JM- 

Proof. Expanding the left-hand side and using (14.11) and the properties of v gives 

we(7, £ ™9) = I n a [2V ( ^ c) n abc ( 7 ) - lbc Yi abc {£ v g)]dvol h . (C.6) 

Jt, 

Applying the Gauss Theorem on the region M~ gives 

ws(7, £»9) = ~ [ V a [2V ( ^ c) n abc ( 7 ) - lbc U abc (£ v g)]dvol g . (C.7) 

J M~ 

The integrand is 

V a (2V ( ^ c) n a6c (7) - lbc Tl abc (£ v g)) = 2V a V (b v c) U abc ( 1 ) 

+ 2V ( ^ c )V a n a6c ( 7 ) - V albc U abc (£ v g) - lbc V a Tl abc (£ v g). (C.8) 
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Using (12.101) and the symmetries of T abcde f , the first and third terms cancel. The remaining 
two terms reduce to 



2V ( ^ c )V a rF 6c ( 7 ) - lhc WJi abc (£ v g) = 2V {b v c) L bc ( 7 ), (C.9) 

where we used the first identity in (j2.12p and L bc (£ v g) = 0. Hence we achieve the desired 
result. □ 

The proof of Theorem IC.ll is completed by combining the results of Lemma IC.2I and 
Lemma El with (JUJ2]). □ 
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